Superconductor-insulator transition in Coulomb disorder 
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Superconductor-insulator transition driven by the decreasing concentration of electrons n is stud- 
ied in the case of the disorder potential created by randomly positioned charged impurities. Electrons 
and Cooper pairs (formed by an non-Coulomb attraction) nonlinearly screen the random potential 
of impurities. Both electrons and Cooper pairs can be delocalized or localized in the resulting self- 
consistent potential. The border separating the superconductor and insulator phases in the plane 
of the concentration of electrons and the length of the Cooper pair is found. For a strong disorder 
the central segment of this border follows the BEC-BCS crossover line defined for a clean sample. 



Superconductor-insulator (SI) transition remains a 
challenging and controversial subject for more than two 
decadesii^i^i^i^i^iii^i^. One way how SI transition with the 
decreasing concentration of electrons n can be envisioned 
is localization of Cooper pairs in a random potential^. 
This approach is good when Cooper pairs weakly over- 
lap and can be considered as repelling each other point- 
like bosons, which at small external disorder experience 
Bose-Einstein condensation (BEC) . At large disorder the 
condensate is fragmented and becomes a Bose insulator. 
If ^ is the size of the Cooper pair created by an (un- 
specified here) attractive interaction between electrons 
the condition of the weak overlap between pairs can be 
written as n^^ ^ 1. 

In the opposite case, n£,^ ^ 1 one may better think 
about separate electrons, which can be localized or de- 
localized by disorder. In delocalized state, the same at- 
traction leads to the Bardeen-Cooper-Schrieffer (BCS) 
superconductivity. Thus, at small enough temperatures 
the metal-insulator transition with the decreasing con- 
centration of electrons leads to the SI transition in the 
fermion picture as well. In this case, however, SI transi- 
tion leads into a Fermi insulator. 

The goal of this paper is to study the zero temperature 
SI transition phase diagram in the plane {S,,n). We sug- 
gest a model in which transition to Bose and Fermi insu- 
lator happens at strongly different concentrations. These 
two segments of the SI border line are connected by the 
long intermediate segment, where n^^ = 1. Here and ev- 
erywhere in this paper we use the scaling approach in the 
large parameter introduced below and drop all numerical 
coefficients. 

The line n^^ = 1 is called BEC-BCS crossover line and 
in a clean sample it is not related to any phase transi- 
tionaiSiii. We, however, show that in our model a long 
segment of BEC-BCS crossover line plays the role of SI 
border line. 

One can not meaningfully describe localization of 
bosons without including their repulsive interaction. 
Therefore, we include the Coulomb repulsion of Cooper 
pairs in our theory. In other words, in the boson limit 
Cooper pairs screen random potential. Similarly, in the 
BCS limit of almost free electrons, the random potential 
is screened by electrons. We assume that the disorder it- 
self is of a the Coulomb origin. Namely, we assume that 



in a three dimensional sample there are randomly dis- 
tributed positive donors with the concentration Nd and 
negative acceptors with the concentration Na, while the 
concentration of electrons is much smaller than both of 
them, n = Nd - Na ^ N, where iV = iV^, + iV^ is the 
total concentration of charged impurities. (This model 
is similar to a heavily doped strongly compensated semi- 
conductor-i^ii^) . The Coulomb nature of both disorder 
potential and interaction between electrons and compact 
pairs makes understanding of screening very easy. 

It is known that on the fermion side of the diagram 
(rt^'^ » 1) even at large concentration of electrons n ^ 
N all electronic states are localized if Na^ <C 1, where 
a — nfi^ /me^ is the effective Bohr radius in the sample, 
K is its dielectric constant and m is the effective mass. 
Thus we assume that the sample is heavily doped. No?" 3> 
1. This of course requires a large k or a small effective 
mass m leading to large a 3> ao, where ao is the lattice 
constant. In this model the disorder is characterized by 
the single parameter Na^ ^ 1. This the parameter which 
is used throughout the paper for scaling estimates. Most 




FIG. 1: The phase diagram of the superconductor-insulator 
transition. On the horizontal axis we plot the length of the 
pair ^ (in units A''"^''^) and on vertical line we plot the elec- 
tron concentration n (in units A'^) both are in logarithmic 
scale. S and I stand for the superconductor and insulator. 
The full line is the border between the superconductor and 
insulator phases. Its sharp corners are results of schematic 
nature of this drawing and actually are rounded. The dashed 
line corresponds to n^"^ — 1. The lower segments of this line is 
crossover between Fermi and Bose insulators, the upper one 
is discussed in the text. 
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of the theoretical models of the SI transition studied in 
literature do not have such a parameter and, therefore, 
have to rely on computer simulations. 

An example, where this model may be applicable, is 
the narrow gap semiconductor PbTe which becomes a su- 
perconductor— with Tc ^ 1.5K, when doping by 1.5% 
of Tl (Pbi_a;Tl3;Te with x ~ 0.015) provides concentra- 
tions of holes up to n ~ 10^" cm""^. The dielectric con- 
stant of PbTe is very large^^, k > 400. As a result the 
Bohr radius a of the Tl acceptor state is so large that 
no freeze out is observed, i. e. for any studied doping 
condition of heavy doping Na^ 3> 1 (where A'' is the con- 
centration of Tl) is fulfilled. It is known that supercon- 
ductivity of PbTe can be induced only by Tl doping. It 
is believedi^ii^ that Tl is essential for attraction between 
electrons which happens via quantum valence fluctua- 
tions of Tl impurities. One, therefore, can imagine using 
Tl concentration to tune attraction and, therefore, the 
pair length ^, while the concentration of holes n can be 
tuned by doping with additional donors, which can com- 
pensate Tl and make n small enough for SI transition. 

Our results are summarized in Fig. [T] The full line 
shows the SI border. It consists of the two horizontal 
lines connected by the central segment of the BEC-BCS 
crossover line n^^ = 1. The lower horizontal line 
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is the border between the superconductor and the Fermi 
insulator phases. There is no dependence on ^ here, be- 
cause electrons are only weakly bound and screen the ran- 
dom potential of charged impurities like free ones. The 
upper horizontal line 
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is the border between the superconductor and Bose insu- 
lator phases. Again ^ is irrelevant, because at such small 
^ the pairs only weakly overlap and can be considered 
as point-like bosons. For a given Na'^ ^ 1 we have got 
n2(A^) ^ ni{N), because bosons have less kinetic en- 
ergy and one needs more bosons for delocalization (see 
derivations of ni{N) and n2{N) below). 

As we mentioned above between the two horizontal 
lines the SI border follows the BEC-BCS crossover line 
n^'^ = 1. Remarkably this part of the border is disorder 
independent. It exists, however, only for a heavy enough 
doping, when Na^ ^1. At Na^ ~ 1 we would get 
n2{N) ~ ni{N) and the intermediate range, where the 
border follows n^^ = 1, shrinks to zero. We see again 
how the whole scaling picture rests on the existence of 
the large parameter Na^ ^ 1. 

The origin of the disorder independent intermediate 
part of the SI border can be understood in the following 
way. Suppose, we cross this border from below, where 
we deal with compact Cooper pairs, which are localized 
by disorder because n < n2{N). Thus, crossing leads 



to almost free electrons with the concentration n larger 
than ni{N). Electrons with such a concentration are 
delocalized. Thus, the crossing of the line n^^ — 1 leads 
from the Bose insulator to the superconductor phase. We 
see that in the presence of a strong disorder the BEC- 
BCS crossover line gets the new meaning. 

Let us make the argument for the intermediate seg- 
ment of the SI border more formal. We know that the 
asymptotic segments n = ni and n — n2 terminate at 
the line n^^ = 1 (the two corners of the border line in 
Fig. [T|). It is natural to assume that between the two 
termination points there is only one intermediate physi- 
cal regime or, in other words, these two points are con- 
nected by a power law n(^) dependence. The only power 
law connecting these points, is, of course, = 1, i.e. 
the line of BCS-BEC crossover. Thus, the SI transition 
border should stick to this line at 112 ^ ri 3> ?t.i. 

Let us now derive the concentrations ni{N) and 
712 (A^). Actually, the concentration ni{N) was derived in 
Ref as the metal-insulator transition in a heavily doped 
strongly compensated semiconductor. It was shown also 
to be in a good agreement with the experimental data 
for compensated semiconductors^"^-. Below we repeat the 
derivation of 77i(Af), because it is a necessary step for our 
derivation of 712 ( A^) . 

Let us divide the sample in cubes with the edge length 
equal R. Due to spacial fluctuations of the concentra- 
tions of donors and acceptor each cube has a random sign 
charge with the absolute value of the order of e{N R^Y^"^ . 
Such randomly alternating charges create the random po- 
tential energy relief of the amplitude 



eV{R) 



kR 



(3) 



This energy diverges at large i?, so that screening even 
by a small concentration of electrons 71 is crucial. To 
discuss this screening let us estimate the characteristic 
fluctuating density SN{R) of charge for fluctuations with 
the characteristic scale R. Clearly 
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The concentration n of electrons can be redistributed be- 
tween wells and hills of the random potential. This re- 
distribution screens all the scales R for which SN{R) = 
{N/R^y/^ < n or, in other words, for R > Rg, where 



Rs = 



1/3 



(5) 



is the nonlinear screening radiu o^^i^^ . All the scales with 
R < Rs remain unscreen, because even when all electrons 
are transferred from all the hills of the potential energy 
to all its wells they are not able to level off the charge 
density of such fluctuations. Since 
thus grows with i?, among remaining scales the most im- 
portant contribution to the random potential is given by 
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R = Rs- Thus, the amphtude of the nonHnearly screened 
random potential energy is 



eV{Rs) 



■T173" 



(6) 



So far we dealt only with the electrostatic energy of elec- 
trons and neglected their kinetic energy. We are talking 
about the limit of zero temperature, so that all the kinetic 
energy is of the quantum origin. Now we should find con- 
ditions when the quantum kinetic energy is small enough 
so that described above regime of localized electrons is 
valid. Clearly the potential energy Eq. ^ is able to lo- 
calize electrons with concentration n if it is larger than 
the Fermi energy of electrons epin) = Ti'r?!^ jlra in its 
wells. In the opposite case e_F(n) eV{Rs) the Fermi 
sea covers the typical maxima of the potential energy and 
the semiconductor behaves as a good metal. Equating 
eV{Rs) and epin) we arrive at the critical concentration 
ni{N) of the SI transition given by Eq. (IT]>i^ii^. 

In the metallic phase electron screening becomes linear, 
the screening radius is given by the standard Thomas- 
Fermi expression Rtf — a/(na^Y^^ and the amplitude 
of the screened potential relief is equal eV{RTF) — 
e^{NRTFy^^ / These quantities match Eqs. (O and 
(|6]) ai n — ni. 

Let us now switch to calculation of the critical con- 
centration 712 (iV) of electrons in charge 2e bosons, corre- 
sponding to the SI transition to the Bose insulator. For 
the insulating phase of localized composite bosons we 
can exactly repeat the above calculation of the nonlin- 
ear screening radius Rs and the random potential energy 
created by screened charged impurities 

The difference between the gas of composite bosons 
and that of fermions becomes important only at the last 
step of calculation of the critical concentration, where 
one has to consider limitations caused by the quantum 
kinetic energy. Many composite bosons can occupy one 
localized level of each potential well. Therefore, condition 
of delocalization of composite bosons is much stronger 
than the condition eV{Rs) < e_F(^^) used for fermions. 
Namely, for delocalization of composite bosons we should 
require that a typical well of the random potential does 
not have a level, or eV{Rs) < h^/mR^. Solving the 
equation 



m 



i?2 



(7) 



for n and using Eqs. ^ and ^ we get the critical con- 
centration of the SI transition Eq. This derivation 
clearly shows why n2{N) ^ ni{N). The potential energy 
amplitude eV{Rs) according to Eq. ^ decreases with 
increasing n. In order to achieve delocalization in the 
composite boson case, we had to make eV{Rs) smaller 
than in the fermion case and this requires the larger con- 
centration 71-2 (iV) 3> ni{N). 

One can arrive at the critical concentration n2 from the 
large n superconductor phase of delocalized Cooper pairs. 



as well. For this purpose we have to start from screening 
in the case when a composite boson can not be localized 
inside the well of the size R and in the first approxima- 
tion its wave function ipQ is the plane wave with a very 
small wave-vector. In the random potential energy with 
characteristic length R and the amplitude 2eV{R) <C 
h^/mR^ this wave function is slightly modulated with 
the small amplitude dip = ■4'o[eV{R) / {fi^ /mR"^)]. For the 
amplitude of the fluctuating density of charge of electrons 
with the scale R this gives 



dn{R) = n 



zV{R) 



{h^/mR^) 



(8) 



All scales R of charge fluctuations for which 5n{R) > 
SN{R) can be screened. Therefore, the equation 5n{R) — 
SN{R) defines the radius Rd of linear screening by delo- 
calized composite bosons. Using Eqs. ([U and (|4]) we get 



Rd 



1/4 



(9) 



(We could not find this result in the literature.) The 
meaning of this screening radius is that the Coulomb 
potential of all fluctuations of charge with R > Rd are 
screened by small changes of wave functions, while all 
the fluctuation scales with R < Rd are unscreened. Us- 
ing Eq. ^ it is easy to see that at 71 ^ ri2 the remain- 
ing amplitude of potential energy fluctuations eV{Rd) ^ 
/mR^, so that we indeed deal with delocalized com- 
posite bosons. Note that as one should expect the lin- 
ear screening radius Rd of the delocalized phase matches 
the nonlinear screening radius Rg of the localized one 
at n = 7i2. Of course, simultaneously V{Rd) matches 
V{Rs). 

One can also verify that on the SI border line in- 
equality Rs ^ which we need to consider screening 
by composite bosons as point like objects fails only at 
^ ~ 7V"1/3(-7Yq3-)2/9^ rpj^jg points resides to the right of 

lower-right corner of SI border line of Fig. [1] This justi- 
fies screening calculations presented above. 

So far we were concerned with the shape of the full 
line separating the insulator and superconductor phases 
on Fig. [T] One may ask whether the two segments of 
the dashed BEC-BCS crossover line n^,^ = 1 abandoned 
by the full line make sense in the insulator phase (the 
right low corner of the diagram) and in the supercon- 
ductor phase (the left upper corner). The answer is that 
these segments of BEC-BCS crossover line imly impor- 
tant crossovers within both phases. 

The dashed line in the insulator phase signals the 
change of a hopping excitation with increasing n from 
a double charged composite boson to an electron. On 
both sides of the dashed line the low temperature variable 
range hopping conductivity obeys the Efros-Shklovskii 
lawil 



a — ao exp 



1/2- 



(10) 



4 



but with two different characteristic temperatures Tq. 
For the large n phase of locahzed individual electrons 
(Fermi insulator) 



where / is the electron localization length and C 2.7 . 
For the low n phase of localized Cooper pairs (Bose in- 
sulator) the double charge 2e replaces e in Eq. pT|) . The 
localization length of a Cooper pair is also smaller than I. 
Thus, the characteristic temperature Tq substantially de- 
creases with increasing n or ^ at the BEC-BCS crossover. 
This leads to the steep increase of the hopping conduc- 
tivity at the dashed line. This phenomenon is similar 
to the reduction of the characteristic temperature Tq by 
the magnetic field in a granular superconductor discussed 
recently in Ref.^^. 

On the other hand, the dashed line segment in the left 
upper corner of our phase diagram signals a change in 
the concentration dependence of the superconductor gap 
and the critical temperature of the superconductor-metal 
transition. At this line both quantities start decreasing 
sharply with the growing n^^. 

Above we considered n and ^ as the two independent 
variables. Strictly speaking this is correct only for n£^^ < 
1. Therefore, small ^ (compact bosons) and intermediate 
segments of the SI border line are found correctly (in the 
scaling sense) . The large ^ segment of our border is above 
the BEC-BCS crossover line n£,^ — 1 and, therefore, ^ 
may be a function of n. (In BCS theory ^ grows with 
n). This, however, can not affect validity of our result 
Eq. ID), because the large ^ segment of the border line is 
not sensitive to ^ at all. This justifies our results for the 
whole SI border. 

Until now we dealt with systems where both impurities 
and electrons reside in three dimensions (we can call it 
3D-3D case). We do not know a 3D-3D system, where n 
can be made tunable in one sample. More promising can 
be the 3D-2D case, where impurities residing in three di- 
mensional space (with the three-dimensional concentra- 
tion N) surround a two-dimensional system of electrons 
(with the two-dimensional concentration n) and affect it 
by the Coulomb random potential. This can be a very 
thin superconductor film on the top of strongly compen- 
sated insulator such as a low quality SiO substrate. In 
such geometry concentration n can be also regulated by 
a gate^^. 

For 3D-2D case one can follow the logic of the calcu- 
lations above and arrive at a similar phase diagram in 
the plane {£,,n). Let us first show the results for this 
case and then outline the calculation. For the critical 
concentrations where localization of electrons and com- 
posite bosons takes place ni and rt2 we arrived at 



where a is the Bohr radius of two dimensional electrons. 
The two horizontal segments of the border line n = ni 
and n = n2 are connected by the BCS-BEC crossover 
line = 1, so that the phase diagram in the plane (n, 
^) looks similar to Fig. [TJ 

The main steps of this calculation are as follows. We 
again start from from cutting three-dimensional space 
filled by impurities in cubes with the edge length equal 
R in such a way that on each side of the plane of the two- 
dimensional electron gas the first layer of cubes touches 
the plane. These two first layers create the "projected to 
the plane" two-dimensional fluctuating charge density of 
impurities {NR^Y^^ /R"^ . Equating this density to n we 
find the nonlinear screening radius R^ = N/n^. Then we 
calculate the amplitude of the nonlinearly screened ran- 
dom potential energy substituting R = Rg into Eq. p|. 
This gives 

N 

eV{Rs) = . (13) 



Equating this energy to the Fermi energy h^n/m yields 
ni, while equating it to h^/mR^ gives n2. In derivation 
of Eq. ([TH)) we followed Refi^ where similar problem was 
studied in order to calculate thermodynamic density of 
states at a small filling factor of a topmost Landau level 
in quantum Hall effect regime in a MOSFET or a GaAlAs 
heterostructure without a spacer. 

In conclusion, we suggested a model with a strong 
Coulomb disorder and strong Coulomb interaction, which 
lets us construct the scaling phase diagram of the SI tran- 
sition in the plane of the concentration of electrons n and 
the characteristic length of the Cooper pair ^ both in 3D 
and 2D systems. The SI border consists of three seg- 
ments. Two asymptotic segments for a large and small 
lengths of Cooper pairs describe transitions in the limits 
of almost free fermions and of compact bosons. They cor- 
respond to the two different concentrations of electrons 
ni and n2. In the strong disorder the threshold concen- 
tration for compact Cooper pairs, 712, is much larger than 
that for almost free electrons, ni. As a result for a strong 
disorder the two asymptotic segments of the border line 
are connected by the third long segment, which follows 
the line of BEC-BCS crossover conventionally defined for 
a clean sampleiSiii. Thus, our model demonstrates that 
the strong disorder can bring about an important impli- 
cation for the BEC-BCS crossover line: it becomes the 
border line of the SI transition. 
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